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Abstract. We generalize the notions of harmonic conjugate functions and Hilbert 
transforms to higher dimensional cuclidean spaces, in the setting of differential 
forms and the Hodge-Dirac system. These harmonic conjugates are in general far 
from being unique, but under suitable boundary conditions we prove existence and 
uniqueness of conjugates. The proof also yields invertibility results for a new class 
of generalized double layer potential operators on Lipschitz surfaces and bound- 
edness of related Hilbert transforms. 
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1. Introduction 

This paper considers higher dimensional analogues of the concept of harmonic 
conjugate functions in the plane. We first review the situation for plane domains. 
Let D be a simply connected domain in R 2 = C. Then, given a harmonic function 
u(z) in D, a harmonic conjugate to u(z) is a second harmonic function v(z) in D 
such that f(z) = u(z)+iv(z) is an analytic function, i.e. satisfies Cauchy-Riemann's 
equations. The function v(z) exists and is unique modulo constants. For example, 
if one requires that v vanishes at some fixed point in D, then we get a well defined 
map u i — > v. Since harmonic functions are in one-to-one correspondence with their 
boundary values, this defines the Hilbert transform 

H D ■ u\ 9D h-> v\q d 

for the domain. The Hilbert transform for a domain D C R 2 concerns only functions 
in D, mapping the real part of an analytic function to its imaginary part. In contrast, 
the Cauchy integral concerns the relation between analytic functions in D and C\D. 
This is best explained through Hardy spaces/projections. Given a function h : dD — > 
C, we form the Cauchy integral 

F(z) = C7 ± /i(z) := / ^-dw, z G D ± , 



2m J 9D w 

where D + := D and D~ := C\D. Taking traces, one obtains two boundary functions 
/ ± (C) := \im z ^^ D ± F(z) such that f + + f~ = hon dD. The Cauchy integral acts 
by projection onto the two complementary Hardy subspaces of the boundary function 
space, consisting of traces / + and /~ respectively. Discarding the exterior Hardy 
function /~ and considering / := / + = u+iv, well posedness of the classical Hilbert 
boundary value problem (BVP) for analytic functions shows that / is in one-to-one 
correspondence with its real part u, as well as with its imaginary part v. (In this 
introduction, we neglect technical details like regularity assumptions on dD and h, 
as well as the fact that maps normally are Fredholm, not exact isomorphisms, in 
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order to explain the main ideas.) From the Hilbert BVP it can be shown that if h 
is required to be real valued, then h is also in one-to-one correspondence with the 
trace of its interior Cauchy integral / = C + h\g D — \h+ 2^P- v -iaD ^( w )/ ( w ~~ z )dw. 
Thus in the diagram 




v 



u 



the double layer potential operator 

Th{z) := Re(C+h\ dD (z)) = \h{z) + -i-p.v. f Im f-^-)h(w), z e dD, 

L Z7T Jq d \W — ZJ 

is an isomorphism. This allows calculation of the Hilbert transform as 
(1) v = H D u = lm(C + (T- 1 u)\ 9D ). 

Note that even though the Cauchy integral, for all D, uses the restriction of l/(w— z) 
as kernel, the kernel of the Hilbert transform depends heavily on D due to the 
factor T _1 . Unfortunately, in the literature the Hilbert transform is often incorrectly 
identified with the Cauchy integral, since it happens to coincide with the (imaginary 
part of) the Cauchy integral for half planes and disks. However, we emphasize that 
for all other domains, these two operators are not the same. 

The aim of this paper is to show the existence and L 2 -boundedness of higher di- 
mensional Hilbert transforms for Lipschitz domains in R n . These Hilbert transforms 
are derived from a Cauchy integral as in ([T]). In the plane, the double layer potential 
operator T is the compression of the Cauchy integral C + \qd to the subspace of real 
valued functions, i.e. T is C + \qd restricted to the subspace of real functions, followed 
by projection back onto this subspace. In higher dimensions, there is a canonical 
Cauchy integral as well, but there are various natural subspaces to compress it to, 
to obtain generalizations of the double layer potential operator. It is of importance 
to establish boundedness and invertibility of such operators, as seen for example in 
(Tjrj. In connection with boundary value problems for Dirac operators, see the works 
|H [21 E] by the first author and Remark EOI for more details, one type of compressed 
Cauchy integrals plays a central role. The compressed Cauchy integrals used in 
this paper to calculate higher dimensional Hilbert transforms have not been studied 
before, to the authors knowledge. 

We next explain the higher dimensional concepts of harmonic conjugate functions 
and Hilbert transforms that we consider in this paper. The plane domain D is 
replaced by a domain in R n and the Cauchy-Riemann system is generalized to the 
Hodge-Dirac system (d + 5)F(x) = 0. Functions F : D — > AR n now in general take 
values in the full 2 n -dimensional exterior algebra 

AR n = A°R n © A x R n © A 2 R n © ... © A n_1 R n © A n R n , 



HILBERT TRANSFORMS AND THE CAUCHY INTEGRAL IN EUCLIDEAN SPACE 



3 



and d and 8 denote the exterior and interior derivative operators 

n 

dF(x) = V a F(x) = tj a djF{x), 

n 

6F(x) = V j F(x) = J2ejj djF(x), 

3=1 

where the exterior and interior products are dual in the sense that (x, Uhy) = (u_ix, y) 
for all multivectors u,x,y G AR™. Multivector fields F taking values in the subspace 
A fe R n we refer to as fc-vector fields. The Hodge-Dirac equation entails a coupling 
between the different fc-vector parts of F. If F — Ff., where F^ : D —>■ A fc R n , 
the differential operators map 
(2) 

S d 8 d 8 d 

■ ■ ■ ~ A fe - 2 R n ~ 1 A fc_1 R n Z A fc R n Z A fe+1 R n 1 ~ A fc+2 R™ ~ • • • 

d 8 d 8 d 8 

so that the Hodge-Dirac equation is equivalent to the system of equations dF k _i = 
—SFk+i, < k < n. If F is monogenic, i.e. satisfies the Hodge-Dirac equation, then 
it is harmonic, i.e. satisfies AF — (d + 5) 2 F = (8d + d5)F = 0, or equivalently each 
of F's 2 n scalar component functions is harmonic. Recall that d 2 = 5 2 = 0. 

Fix < k < n and consider a harmonic k- vector field U : D — > A fc R n . We say 
that V\ : D — > A fe_2 R n and V2 : D — > A fc+2 R n form a pair of harmonic conjugates 
to U if (d + <5)(X4 + U + V 2 ) = 0, or equivalently dU = -SV 2 , 5U = -dV 1 and 
dV 2 = SVy = 0. 

Example 1.1. When n = 2 and k = this reduces to the classical situation. Indeed, 
consider an analytic function / = u + iv. We identify R ~ A°R 2 and iR ~ A 2 R 2 . 
Fixing an ON-basis {ei, e 2 } for R 2 , we identify i ~ ei a e2- Then 

(d + 5){u + if) = (ei a d\u + e 2 a <9 2 w) + (ei j zSif + e 2 j i9 2 w) 

= (<9i-u - 9 2 f )ei + (d 2 u + div)e 2 = 

coincides with the Cauchy-Riemann equations. Hence v is a classical harmonic 
conjugate to u if and only if V 2 = v e± a e 2 : D — > A 2 R 2 is a harmonic conjugate to 
u in the sense of the Hodge-Dirac system. (In this case, the harmonic conjugate V\ 
vanishes.) 

In the general case, we observe that a necessary condition for such Vi, V 2 to exist 
is that U is two-sided harmonic, i.e. 5dU = = dSU. We also observe that V±, V 2 
are only well defined modulo two-sided monogenic fields, i.e. the differences V\ — V{ 
and V 2 — V 2 of two sets of harmonic conjugates satisfy <i(Vi — Vf) = = 5(Vi — V(), 
i — 1,2. When 1 < A; < n — 1, the two-sided monogenic /c-vector fields form an 
infinite dimensional space (see Corollary 13.111) . Thus, in order to obtain a uniquely 
defined higher dimensional Hilbert transform, further conditions need to be imposed 
on V\ and V 2 , so that there is a well defined map 

U^Vx,V 2 . 

In this paper we consider one possible such further condition on harmonic conjugate 
functions, which extends the above technique of calculating conjugates with the 
Cauchy integral and double layer potential operators to higher dimension. Under 
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this further condition we say that the harmonic conjugate functions are of Cauchy 
type. Since all component functions of U, V\ and Vi are harmonic, these fields are in 
one-to-one correspondence with their trace on 3D. Thus, equivalently we will have 
a Hilbert type transform U\qd | — ► ^i|aD, V^ao for D. 

The outline of the paper is as follows. In Section [2] we introduce the higher 
dimensional Cauchy integral associated with the Hodge-Dirac system, and prove 
existence and uniqueness results for Cauchy type conjugates to scalar functions, i.e. 
k = or k = n. This amounts to proving invertibility of the classical double layer 
potential operator for domains in R n . Both the boundedness and invertibility of 
this singular integral operator on L p -spaces on Lipschitz boundaries (1 < p < oo for 
boundedness and 2 — e < p < oo for invertibility) are deep results, but are by now 
well known facts. 

Section H] contains the main new results of the paper, Theorem 14. 11 and establishes 
existence and uniqueness of Cauchy type conjugates to k- vector fields, 1 < k < n — 1. 
This general case is more involved than the scalar case, since the generalized double 
layer potential operators which appear will not in general be invertible, not even 
Fredholm, as they have infinite dimensional null spaces and cokernels. However, us- 
ing the theory of boundary value problems for Dirac operators (which is reviewed in 
Section[S]), we manage to show invertibility of the operator acting from a complement 
of the null space to its range, in a natural L 2 -based Hilbert space. 

In the final Section [5] we illustrate the non- uniqueness of harmonic conjugate 
functions in higher dimensional euclidean spaces by constructing different conjugate 
functions which are not in general the Cauchy type conjugates. This second con- 
struction is based on the theory of Hodge decompositions, and the obtained harmonic 
conjugate functions are said to be of Hodge type. 

In the literature, various generalizations of harmonic conjugate functions to higher 
dimensional euclidean spaces can be found. A classical generalization for the upper 
half space, using divergence and curl free vector fields, was introduced in harmonic 
analysis by Stein and Weiss [11], see Stein [12J and Example 12. 7( 1). A generalization 
more similar to our construction is due to Arzanyh pQ, who studied two-forms B : 
D — > A 2 R 3 conjugate to scalar functions in three dimensional space. In the setting 
of Clifford analysis, without dealing directly with the more fundamental differential 
operators d and 5, there is work on Hilbert transforms and harmonic conjugate 
functions in euclidean space by Brackx, De Knock, De Schepper and Eelbode [6]. 
See also the references therein for calculations on special domains like the unit ball. 



Writing D := d+5 for the Hodge-Dirac operator, this is an elliptic first order par- 
tial differential operator whose square is the Hodge-Laplace operator D 2 = A. Just 
like the exterior and interior differentiation use the exterior and interior products, 
the Hodge-Dirac operator uses the Clifford product as 



The Clifford product is the unique associative algebra product a on AR n , with 
identity 1 e A°R n = R, such that 



2. Hilbert transforms for scalar functions 



n 




(3) 



V AW = V JW + V AW, 



W Atl = W Lt) + W AV 
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for all vectors v G A x R n = R n and all multivectors w G AR™. The main difference 
between the complex product in R 2 and its higher dimensional analogue in R n , the 
Clifford product, is that the latter is non- commutative. Here j and l denote the 
left and right interior products, defined as the operations adjoint to left and right 
exterior multiplication, i.e. 

(4) (w j x, y) — (x, w a y), (x i_ w, y) — (x, y a w), w,x,y E AH n . 

Following standard notation, we write w± AW2 ='■ W1W2 for short. Important to this 
paper is the following mapping property of the Clifford product. If v G A 1 R™ = R™ 
and w G A fc R n , then 

vw G A fc - x R n © A fe+1 R n . 

This is clear from 

Concretely, if {ej}™ =1 denotes the standard ON-basis for R™, the induced ON- 
basis for the space of fc-vectors A fc R n is {e s }\ s \=k, and in total {e s } sC n is the induced 
basis for AR n , where n := {1,2,..., n}. For a subset s = {si, . . . , Sk} C n, where 
1 < $1 < ■ ■ ■ < Sk < n, we write e s := e Sl a . . . a e Sk . In this induced basis the 
exterior, interior and Clifford products are 

e s a e t = e(s, t) e sUt , if s Pi t — 0, and otherwise 0, 

e s je t = e(s, t \ s) e t \ s , if s C t, and otherwise 0, 

e s l e f = e(s \ t, t) e s \ t , if t C s, and otherwise 0, 

e s e t = e(s,t) e tAs , 

where the permutation sign is e(s,t) := (— 1)I{( s ''*j) Gsx ' ; s i>*j} = ±1 anc l s A t := 
(s \ t) U (t \ s) denotes the symmetric difference of index sets. 
The radial vector field 

E{x) := —y^-, x G R", 

0"n-l F| n 

where u n -\ is the area of the unit sphere S" --1 , is divergence and curl free for x ^ 0, 
and is a fundamental solution to the Hodge-Dirac operator. Using the associativity 
of the Clifford product, Stokes' theorem gives a Cauchy type integral formula 

(5) F{x) = [ E(y~ x)u(y)F(y)da(y), x G D, 

JdD 

for monogenic fields F : D — > AR n , i.e. solutions to DF = 0. Here v{y) denotes 
the unit normal vector field on 3D pointing outward from D, and da denotes scalar 
surface measure on 3D. Note that u(y) must be placed between the two factors 
because of non-commutativity. 

Example 2.1. If n = 2 and F : R 2 w C -> C w A°R 2 © A 2 R 2 , the Cauchy formula 
([5]) reduces to the standard one. Indeed, 

r-i/ 1 dw 

E(y - x)u(y)da{y) = — 

2-ki w — z 

if we identify % = eie 2 , w = eiy, z = t\X and dw/i = e\vdo~. 

With these algebraic preliminaries, we next turn to analysis. For estimates, we 
use the notation X w Y to mean that there exists a constant C, independent of the 
variables in the estimate, such that C~ X Y < X < CY . Similarly X < Y means that 
X < CY. 
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To avoid topological technicalities, we shall restrict attention to the following two 
types of strongly Lipschitz domains C R n (i.e. domains whose boundaries are 
locally graphs of Lipschitz functions). We write £ := dD + = dD~ for the Lipschitz 
boundary between D + and D~ = R n \ D . By D we denote either D + or D~ . 

• A graph domain D + = {x ; x n > (j){x\, . . . , n n -i)} above the graph of a 
Lipschitz regular function : R n_1 — > R. Here D~ denotes the domain 
below the graph. 

• An interior domain D + , being a bounded domain which is Lipschitz dif- 
feomorphic to the unit ball, and whose boundary is locally the graph of a 
Lipschitz function (in suitable ON-bases). The exterior domain D~ is the 
interior of the unbounded complement of D + . 

The unit normal vector field u(y) on £ is always assumed to point into D~ , i.e. 
the region below the graph or into the exterior domain. We define non-tangential 
approach regions 7(1/) C D, y G E, for these Lipschitz boundaries. For graph 
domains D , fix c\ greater than the Lipschitz constant for E, and let 

7(?/) =7(y,# ± ) := {(x',x n ) G R^ 1 x R ; ±(x n - y n ) > Cl \x'-y'\}, 

for y = (y', y n ) G E. For exterior and interior domains, and y G E, consider the 
coordinate system around y in a neighbourhood of which E is a Lipschitz graph. 
The approach region 7(7/, D^) is defined as the truncated part of the cone, where 
dist (x, y) < c 2 and c 2 denotes a sufficiently small constant. 

The boundary function spaces we use are the spaces £ p (£; A), where A := AR™. 
For a field F in D^^, define its non-tangential maximal function 

N*(F)(y):= sup \F(x)\, y G E. 

A fundamental theorem in harmonic analysis and singular integral theory due to 
Coifman, Mcintosh and Meyer [7], states that the Cauchy integral is bounded on 
L P (E; A) on any Lipschitz surface E. By surface, we shall mean a hypersurface in 
R n . 

Theorem 2.2. Let be Lipschitz graph, interior or exterior domains, and fix 
1 < p < 00. Let h G £ p (£; A) and define the monogenic field 

C+hix) :=± [ E(y - x)u(y)h(y)da(y), x G D ± . 

Then \\N *(C + h)\\ p + 1 1 A^"^ ( C ~ ) 1 1 ^, < C||/i|| p for some C < 00 depending only on p 
and the Lipschitz constants for the graphs describing E. 
The principal value Cauchy integral 

Eh(x) := 2p.v. / E(y — x)u(y)h(y)da(y), x G E, 
is 

exists a.e. and defines a bounded operator E : L p (E; A) — > L P (E; A) such that 
E 2 = I. The boundary traces f + {z) := \ivn-x^z,xe-y(z,D+) C + h(x) and f~~(z) := 
\im x ^ z ,x£'y{z,D-)C~h(x) exist for a. a. z G E and in L p , and 

E + h:= f + = \{h + Eh) and E~ h := /~ = \{h - Eh) 

define L p -bounded projection operators. 
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Let h G L p (Yi] A k ) be a A;-vector field and consider harmonic conjugates in D + . 
The Cauchy integral produces a monogenic field 

F = Vi + U + V 2 = C + h : D + — > A k - 2 R n © A fc R" © A fc+2 R" 

in D + . Indeed, the mapping properties of the Clifford product show that multi- 
plication with the normal vector gives vh : E — > A 1 R n © A fe+1 R n , and a similar 
splitting when multiplying with the vector E(y — x) shows that 

(6) C + h(x) = J E(y — x) j (y(y) j %)) da(y) 

+ ^ - x) a j h{y)) + £(y - x) j a %))) da{y) 

+ J^E(y-x)A v{y) a h{y) da(y) = V 1 + U + V 2 . 

Definition 2.3. Given a two-sided harmonic field U : D ± — > A fc R n , i.e. <5<i?7 = = 
c/5f/, we say that Vi,V2 are Cauchy type harmonic conjugates to Z7 if there exists 
/i : E -> A fc R n such that U = (C7 ± /i) fc , ^ = (C7 ± / i ) fc _ 2 and 7 2 = (C ± /i) fc+2 , where 
subscript denotes the /c-vector part of a multivector. We call h the dipole density 
of the system ]/]_, U, V 2 of harmonic conjugate functions. The corresponding map of 
boundary values 

£/|E^Vi| s ,V 2 | s , 

we refer to as the (Cauchy type) Hilbert transform for the domain D. 

The following theorem on existence and uniqueness of Cauchy type harmonic 
conjugates to scalar functions (k = 0, n) is the main result of this section. 

Theorem 2.4. Let D C R" be a Lipschitz graph, interior or exterior domain and 
assume that 2 < p < oo. 

(i) Let U : D — > R = A°R n be a harmonic function such that N*(U) G L p (E). 
// .D is an exterior domain, also assume that lim x .^ 00 U = and /ias trace 
u = U\s such that J uip da(y) = 0, where tp is the function from Theorem \2.5i 
Then there is a unique Cauchy type harmonic conjugate V = V 2 : D — > A 2 R™ 
to U , and a dipole density h G L P (E) ; swc/i t/iat 

ll^(^)L + ll%<ll^(^)ll P - 

// D is a graph or an interior domain, then h is unique, and if D is an 
exterior domain, then h is unique modulo constants. 

(ii) In the case k = n, (i) remains true when U : D — ► A°R n = R and V = V 2 : 
D -> A 2 R n are replaced by U : £> -> A n R" w R anc? V = Vi : D -> A"~ 2 R" . 

The scalar cases = and k = n are significantly more straightforward than 
the non-scalar case 1 < k < n — 1 (to be treated in Section HJ) as they reduce to 
the question whether the classical double layer potential equations are invertible, as 
explained in the two dimensional case in the introduction. On E, define the principal 
value double layer potential operator 

(7) Kh(x) := 2p.v. J (E(y - x),n(y))h(y) da(y) = (Eh(x)) Q , 

for h : E — » R, x G E. The boundedness of K in L p (E), 1 < p < oo, is a direct 
consequence of Theorem l2.2l Invertibility of I±K on the other hand, which the proof 
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of Theorem [231 uses, is not true for all p on a general Lipschitz surface. Invertibility in 
L 2 (S) was proved by Verchota [13J via the method of Rellich estimates. Invertibility 
in the range 2 < p < oo was proved by Dahlberg and Kenig [8j by atomic estimates 
in real Hardy space duality and interpolation. 

Theorem 2.5. Assume that 2 < p < oo and let E be a Lipschitz graph or the 
boundary of an interior / exterior domain. Then 

I + K : L P (E) — ► L P (E) : /i i — ► h + (Eh) = 2(E + h) 

is an isomorphism. This is also true for I — K : h i— > (E~h)o in the case of a graph 
domain. In the case of an exterior domain, I — K is a Fredholm operator with null 
space consisting of constant functions and range consisting of all u G L P (H) such 
that L uip = for some ip G L q {T), where 1/p + 1/q — 1. 

We remark that Theorems 12.21 12~31 and 12751 can be generalized to strongly Lipschitz 
domains with more complicated topology. In this case existence and uniqueness of 
conjugates hold only modulo finite dimensional subspaces, and I + K are Fredholm 
operators with higher dimensional null spaces and cokernels. 

Proof of Theorem \2.4\ To prove (i), take h G L P (E; A ) and define U := (C ± /i)o and 
u := U\y, = |(J ± K)h G L p (H; A ). Theorem 12.51 determines uniquely h such that 
u — u, possibly modulo constants in the case of an exterior domain. In any case, this 
defines uniquely a Cauchy type harmonic conjugate V := (C ± h) 2 , since C~ maps 
constants to zero in an exterior domain. 

To prove (ii), consider a system of Cauchy type harmonic conjugates C ± h = V+U, 
where h : E —> A n R n , V : D — > A n ~ 2 R n and U : D —> A n R n . Introduce the operator 
U i— > Ue-n (he. the Hodge star operator for differential forms up to a sign) which maps 
k- vectors ton - k- vectors. We have C ± {hen) = Uen + Ven, where the functions take 
values in A°R ri , A°R ri and A 2 R n respectively. This reduces (ii) to (i), since U i— > Ue-n 
commutes with the Dirac operator and the Cauchy integral by the associativity of 
the Clifford product. □ 

We note from the proof the following relations between Cauchy type harmonic 
conjugate functions and Hilbert transforms of scalar functions, the Cauchy integral 
and the classical double layer potential operator. 

Corollary 2.6. Let V : D + — > A 2 R" be the Cauchy type harmonic conjugate to the 
harmonic function U : D + — ► A°R n = R, with suitable estimates of non-tangential 
maximal functions. Then 

U + V = 2C + ((J + K)- l u), 

where C + is the (interior) Cauchy integral, K is the double layer potential operator, 
and u = U\s- Taking the trace v = V\y, of the conjugate function, the Cauchy type 
Hilbert transform of u is 

U \ — > v = (I + E)(I + K)- l u. 

Replacing C + , I + K and I + E with C~ , I — K and I — E, the corresponding 
formulae hold for the domain D~ . 
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Example 2.7. (1) If D = R™ is the upper half space with the flat boundary 
E = R™" 1 , then K = since E(y — x) in this case is orthogonal to v(y). Thus 
h — 2u and 

V(x) = 2(C + u) 2 (x) = e n a — [ JL^- U (y) dy, x e R n + , 

&n-l jRn-l \V - X\ n 

where e n is the basis vector normal to R n_1 , so in this case the Hilbert transform 
U\-Rn-i i — > V|h.»-i for the upper half space coincides with the bivector part of the 
principal value Cauchy integral u i— > {Eu)i- 

A classical higher dimensional notion of harmonic conjugates, using divergence 
and curl free vector fields in the upper half space R™, is due to Stein and Weiss [TT] . 
and we refer to Stein's book [12] for further details. The upper half space has the 
special property that the vector e n normal to <9R™ = R n_1 is constant. Split a 
vector field F in R™ into normal and tangential parts as 

F(x) = U{x)e n + U(x), 

where U is a scalar function. Stein and Weiss consider the tangential vector field 
U as a harmonic conjugate to U, if F is a divergence and curl free vector field, i.e. 
if F is monogenic. Since the Clifford product is associative, this is equivalent to 
U + Ue n : R™ — > A°R n © A 2 R" being monogenic since 

D(C7 + Ue n ) = B((Ue n + U)e n ) = (D(Ue n + U))e n = 0. 

Due to the very special geometry of R™, the bivector field Ue n will in fact be 
the Cauchy type harmonic conjugate to U. Indeed, if V denotes the Cauchy type 
conjugate, then as noted in the introduction, the difference Ue n — V is a two-sided 
monogenic bivector field. Moreover V(x) = e n A J Rn E(y—x)u(y)dy, since the normal 
vector is constant. Hence e n a (Ue n — V)|rt.-i = 0. From Theorem 13.31 below we 
deduce that Ue n — V — 0, since there are no non-trivial monogenic field which are 
normal (in the sense of Definition 13. lft on the boundary. 

Thus, from the above relation Ue n = 2(C + m)2, the Stein- Weiss tangential vector 
field U is seen to be 

U(x) = — ^— / \-u(y)dy, x = (x',x n ) G R™, 

o-n-i Jr« \x - y\ n 

and taking the trace of U, the n — 1 component functions of U\nn-i are the Riesz 
transforms of u. We remark that harmonic conjugates in the sense of Stein and 
Weiss do not generalize to more general domains D, since they depend on a canonical 
direction e n . 

(2) If D is the unit disk in the plane, then 

Kh(x) = -p.v. / ( f ~ X 'l } h(y) da{y) = [h], 

is \y-x\ 

where [h] denotes the mean value of h, regarded as a constant function. This gives 
(J + K)~ l u = u - [u]/2 and 



V(x) = (2C+(J + K)- l u) 2 (x) = (2C + u - [u]) 2 = 2(C + u) 2 {x) 

1 f 27T 

-u(y) da{y) — ► (Eu) 2 {e lt ) = ^P-V-J cot((t - s)/2)u(s) ds, 



Js \y - x\ 
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when x — > e lt , y = e ls , i = eie 2 , and the limit is pointwise a.e. and in L p , 1 < p < oo. 
Also here the Hilbert transform coincides with the imaginary part of the principal 
value Cauchy integral. Note that this is far from being the case for more general 
domains, not even when n = 2, or for higher dimensional spheres as we shall see 
below. For the exterior of the unit circle, we see that (I — K)h = h—[h]. Thus I — K 
has constants as null space and its range consists of all function with mean value 
zero. Hence the function ip in Theorem 12.41 orthogonal to the range, is a constant 
function. 

(3) For the unit sphere E in R n , K is the operator 

1 f h(y) 

This is weakly singular and therefore compact on all L p spaces. In fact, this is 
true whenever E is a smooth surface, since the normal vector will be approximately 
orthogonal to E(y — x) in the kernel when y is close to x. Through a limiting 
argument, it was proved by Fabes, Jodeit and Riviere [9] that K is a compact 
operator on L p , 1 < p < oo, whenever E is a bounded C 1 regular surface. 

3. DlRAC BOUNDARY VALUE PROBLEMS 

In this section we describe parts of the operator theoretic framework for Dirac 
boundary value problems, developed by the first author in his PhD thesis [3], which 
we use in Section H] to extend Theorem 12. 41 to general k- vector fields. In this section, 
we focus on explaining the main ideas of proofs, but give references to full proofs. 
By N(T), R(T) and D(T) we denote the null space, range and domain of an operator 
T. 

The basic picture is that the boundary function space 

(8) L 2 (E) := L 2 (E; A) = L 2 (E; A ) © L 2 (E; A 1 ) © ... © L 2 (E; A"" 1 ) © L 2 (E; A n ) 

splits in two different ways into pairs of complementary closed subspaces 

L 2 (E) = E + L 2 © E-L 2 = N + L 2 © N~ L 2 . 

In the first splitting L 2 = E + L 2 © E~L 2 , the subspaces E ± L 2 denote the Hardy 
type subspaces associated with the Dirac equation, i.e. E + L 2 consists of traces 
of monogenic fields in D + and E~L 2 consists of traces of monogenic fields in D~ 
which vanish at infinity. The Hardy subspaces E ± L 2 = R(-E ± ) are also the ranges 
of the Hardy projection operators E ± in L 2 (S) from Theorem 12.21 which explains 
the notation. There is a one-to-one correspondence between / = F\-£ G E ± L 2 and 
F = C ± f : —>■ AR n , and we sometimes identify F and /, referring to F as 
belonging to the Hardy type subspace. 

In the second splitting L 2 = N + L 2 © N~L 2 , which is pointwise, the subspace 
N + L 2 consists of all fields tangential to S, and N~L 2 consists of all fields normal to 
S. 

Definition 3.1. A multivector field / : E — > AR n is tangential if v{x) j f(x) = 
for almost all x G E, and it is normal if u(x) a f(x) = for almost all x G E. 

The two projection operators N ± are 

N + g := v j {y a g) and iV~ g := v a [y j g) . 
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This tangential/normal splitting is orthogonal, so that \\g + + (?-|| 2 = ||<7+|| 2 + ||<7-|| 2 , 
when g± G N ± L 2 . On the other hand, the Hardy space splitting is not orthogonal, 
only topological in the sense that ||/ + + /_|| ~ ||/ + || + ||/_|| when f± G E ± L 2 . 

The operator theoretic problem underlying boundary value problems is to under- 
stand the relation between the splitting E + L 2 © E~L 2 and the splitting N + L 2 © 
N~L 2 . 

Example 3.2. Consider the following basic Dirac BVP consisting in finding F : 
D + — ► AR n solving the Dirac equation (d + 5)F = in D + with given normal part g 
on the boundary E. Under appropriate regularity assumptions, this means exactly 
that we are looking for / = F\^ G E + L 2 such that N~ f = g. Uniqueness and 
existence of such /, for each g, is clearly equivalent to the restricted projection 

N~ : E + L 2 — ► N~L 2 

being an isomorphism. 

In general, there are topological obstructions preventing N~\e+l 2 from being an 
isomorphism. However, modulo finite dimensional spaces, the operator is always 
invertible. 

Theorem 3.3. Let E be any strongly Lipschitz surface. Then the restricted projec- 
tion N~ : E + L 2 — > N~L 2 is a Fredholm operator of index 0, i.e. has closed range 
and finite dimensional kernel and cokernel of equal dimensions. The same is true 
for all eight restricted projections 

N + : E ± L 2 — ► N + L 2 , N~ : E ± L 2 — ► N~L 2 , 

E + : N ± L 2 — > E + L 2l E~ : A^ ± L 2 — > E'L 2 . 

If H is a Lipschitz graph, then all these maps are isomorphisms. 

The key ingredient in the proof is a Rellich type estimate. The strong Lipschitz 
condition on S, i.e. that S is locally the graph of a Lipschitz function, shows the 
existence of a smooth vector field 9 which is transversal to S, i.e. (v(x), 9(x)) > c > 
for all x G E. Basic identities for the Clifford and interior products show that 

|/|>, 9) = \f\ 2 \(v9 + 9v) = (fv, fO) = (-2u jf + vf, fO), 

where / i— > f is the automorphism which negates fc-vector fields with odd k. Thus 
an application of Stokes' theorem yields the following Rellich identity 

(9) / \f\ 2 (u,9)da(y) = -2 [ (v -> fJO) da(v) + Y. [[ (F,ejF(d s 6)) dx, 

JT, JT, ■ 1 J JD+ 

for all / = F|s G E + L 2 , and therefore the estimate ||/|| < ||A^ _ /|| + HFH^fsuppe)- If 
E is a Lipschitz graph, we can choose 9 = —e n , in which case the last term vanishes 
and it follows that the restricted projection iV~ : E + L 2 — > N~L 2 is injective and 
has closed range. More generally, the map / i— > F in the last term in the estimate 
can be shown to be compact, from which it follows that iV~ : E + L 2 —> N~L 2 has 
finite dimensional null space and closed range. Finally the index of the restricted 
projection can be shown to be zero through either a duality argument or the method 
of continuity. For details we refer to [2]. 

As we shall make frequent reference to it, let us state the well known method of 
continuity. 
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Theorem 3.4 (Method of continuity). Let X and y be Banach spaces, and assume 
that T\ : X — > y, A G [0, 1] ; is a family of bounded operators depending continuously 
on A. IfT\ are all semi-Fredholm operators, i.e. has closed range and finite dimen- 
sional null space, then the index, i.e. dim(y/R(T\)) — dim N(T\), of all operators 
T\ are equal. 

Another operator in L 2 (E) of importance to us, besides E^ 1 and N^, is the fol- 
lowing unbounded first order differential operator T. 

Definition 3.5. Let E be a strongly Lipschitz surface. Denote by T the unique 
closed operator in L 2 (E) with dense domain D(r) C L 2 (E) and the following action. 

(i) If / = /+ + /-= F + | s + F_| s G E + L 2 © E~L 2 in the Hardy space splitting, 
then T acts by exterior (= — interior) differentiation on the monogenic fields 
F ± in D ± as 

17 = (dF + )| E + (dF_)| s = (-<fF + )| s + {-5F_)\v. 

(ii) If / = fx + v a f 2 G N + L 2 © N~L 2 in the tangential/normal splitting, where 
A) /2 € N + L>2, then T acts by tangential exterior and interior differentiation 
on the two parts respectively as 

rf = dvf 1 + VA(8vf 2 ), 

where ds and 5s denote the intrinsic tangential exterior and interior differ- 
entiation operators on the surface E. 

Recall that a bilipschitz parametrization p : R rt_1 — > S, locally around a point 
y G E, induces a pullback p*, mapping tangential multivector fields N + L 2 to 
L 2 (R n_1 ; AR™ -1 ). Exterior differentiation commutes with this pullback, i.e. ds/ = 
(p*) _1 d R n-ip*/. Dual to this, a reduced pushforward p* in the terminology of [3J, 
intertwines 5s and 5r«-i. 

For Definition 13.51 to make sense, one needs to show that the operators in (i) and 
(ii) coincide. This is a consequence of the following proposition. 

Proposition 3.6. The following intertwining relations for exterior and interior dif- 
ferentiation operators hold. 

(i) IfU : D ± -> AR n and N*(U),N*(dU) G L a (E), tfien iV+[/|s G D(d E ) and 
r(A^+C/|s) = d E (iV + t/|s) = iV+(dC/| E ). If U : D ± ^ AR n and /V*([/) ; 
iV,(5C/) G L 2 (E) ; i/ien f j £/|s G D(5 E ) and r(JV-tf| E ) = v a 5 s (i/ j 17|e) = 
iV-(5C/| s ). 

(ii) If h £ N + L 2 , d-zh G L 2 (E) and x ^ E, £nen 

d / E(y - x)v{y)h{y) da{y) = / E(y - x)v(y)(d s h)(y) da(y). 



If h G N L 2 , 5s(^ j /i) G £ 2 (E) and x ^ E ; iaen 




The trace result iV + (d?7|s) = ds(iV + [/|s) in (i) is a special case of the fundamental 
fact that the exterior differentiation and pullbacks commute. Indeed, if i : E — > R" 
denotes inclusion, then iV + [/|s = i*(U). The trace result iV~(5[/|s) = ^ a 5 E (i> j 
C/|s) can then be obtained by Hodge star duality. For more details of the proof of 
Proposition 13.61 we refer to jH Proposition 4.10] and [3J Proposition 6.2.5]. 
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The usefulness and relevance of the operator T to this paper, is that it is closely 
related to the decomposition (jSJ), as the following proposition shows. 

Proposition 3.7. The operator V relates to the splitting o/AR" into homogeneous 
k-vectors as follows. 

(i) If he L 2 (E; A fc ) and T(N + h) = 0, then (c7 ± /i) fc+2 = 0. If he L 2 (E; A k ) and 
r{N~h) = 0, then (C7 ± /i) fc _ 2 = 0. Thus, if h e L 2 (E; A fc ) and Th = 0, then 
C ± h is a k -vector field. 

(ii) For a monogenic field F = ^^ =0 -Ffc in D ± , i.e. (d + 8)F = 0, where Fk : 

— >• A fc R n , the following are equivalent. 

— All homogeneous parts Ft of F are monogenic, i.e. (d + 5)Fk = 0. 

— F is two-sided monogenic, i.e. DF = = FT), where FD := Y^j=i{®jF) e j- 

— F satisfies dF = 5F = 0. 

The result (i) shows that T governs the off-diagonal mapping of the Cauchy in- 
tegral, i.e. the first and last terms in (jSJ). Indeed, an integration by parts rewrites 
these terms as single layer potentials 

E(y - x) Au(y) hh{y)da{y) = / $(y - x) v{y) a (Th(y)) da(y), x^T,, 
s Jt, 



E(y - x) j (y(y) j h(y)) da(y) = - / $(y - x) v{y) j (Th(y)) da(y), x £ E, 

where $>(x) denotes the fundamental solution for the Laplace operator and E(x) = 
V$(x). For proof of Proposition 13.71 we refer to [21 Lemma 4.13, Proposition 4.5]. 

Since F acts by exterior and interior differentiation, it is clear that T 2 = 0, or 
more precisely R(r) C N(r). We have inclusions of function spaces 

L* (£) C ^(S) c L 2 D (S) c L 2 (E), 

where Lf (E) := R(r), Lf(E) := N(r) and Lf (E) := D(r). Here Lf(E) is always a 
closed subspace of L 2 (E) and L 2 :) (E) is always a Hilbert space densely embedded in 
L 2 (E). The domain -^(E) is equipped with the graph norm ||/|||, := ||/|| 2 + ||r/|||, 
which makes it a Hilbert space. The range Lf("E) is equipped with the range norm 



^:=inf{||r M || 2 +||n|| 2 ; u G D(r), Yu = /}, 

Rl 



which makes it a Hilbert space. The properties of the range L 2 (E) depends on the 
surface E, as the following lemma shows. 

Lemma 3.8. If E is an unbounded Lipschitz graph, then L R (Y1) is dense and not 
closed in (E). If T, is a bounded Lipschitz surface, then L R (E) is a closed subspace 
of L>2 (S) of finite codimension. In particular, if D + is Lipschitz diffeomorphic to 
the unit ball, then the codimension is 4 and 

L R (X) = {/ e Lf(E) ; /o = / n = / UAf n ^da= [ u^f 1 d<T = 0}, 

where : E — >• A fe R n denotes the k-vector part of f . 

Proof. Since N + L 2 and N~L 2 are invariant under T, we may consider tangential and 
normal multivector fields separately. Moreover, the two operators d and 5 acting in 
R" satisfy 

5(Fen) = (dF)e w , 
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and taking the normal part of the trace of this identity, Proposition 13. 6( i) shows that 
r(/en) — (r/)e„ for all tangential / G N + L 2 . Thus the actions of T on N + L 2 and 
N~L 2 are similar, and it suffices to consider T = acting on tangential multivector 
field. Here the stated results are well known facts from de Rham cohomology. □ 

Instead of working with the projections E ± and A r± , it is often convenient to work 
with the associated reflection operators E := E + — E~ and iV := iV + — N~, where 
E 2 = N 2 = I, which explains the naming. Here E is the principal value Cauchy 
singular integral from Theorem 12.21 and 

Nf = vfv 

is the operator which reflects a multivector field across E. Following the boundary 
equation method developed in [31 Hj, we shall make use of the rotation operator 

EN f{x) = p.v. / E(y - x)f(y)u(y)da(y). 

The important connection between EN and the restricted projections above, is that 

I + EN = 2(E + N + + E~N-), N(I + EN)N = 2(N + E + + N'BT), 

I - EN = 2(E + N~ + E~N + ), N(I - EN)N = 2(N + E~ + N~E + ). 

For example, this shows that / + EN is the direct sum of the restricted projections 
E + : N + L 2 — > E + L 2 and E~ : N~L 2 — > E~L 2 . Thus, in order to prove that all 
eight restricted projections are Fredholm operators, it suffices to prove that the two 
operators / ± EN are Fredholm operators on the full space ^(E; A). We record 
the following generalization of Theorem 13.31 which was proved in [2, Theorem 4.15] 
through Rellich estimates involving a pair of monogenic fields F^ 1 : — > AR n and 
the method of continuity. 

Theorem 3.9. Let E be a strongly Lipschitz surface. Then X + EN : L 2 (E) — ► L 2 (E) 
is a Fredholm operator with index zero for all X G R (and more generally in a double 
sector around the real axis). 

The last result we shall need is the following analogue of Theorem 13.31 for the 

subspaces Lf(E), L 2 (E) and Lf(E). 

Theorem 3.10. All four projections E^ and N ± leave each of the subspaces L 2 (E), 
L 2 (E) and L 2 (Ju) invariant and act boundedly in them. All eight restricted projec- 
tions 

N + : E ± L 2 — ► N + L x 2l N~ : E ± L X 2 — ► N~Lf, 

E + : N^L^ — > E + L x 2 , E~ : A^L* — > E~ L x 2 

are Fredholm operators, for x = R, N, D . All eight maps are injective when x = R, 
i.e. when acting in the range L^(E), for all strongly Lipschitz surfaces E. 

Proof of Theorem \3.1(A That E ± and N^ 1 act boundedly in all three subspaces is 
clear from (i) and (ii) in Definition 13.51 As noted above, the Fredholm property of 
all eight restricted projections will follow if we prove that / ± EN are Fredholm 
operators on L|(E). 

(1) Fredholmness of the operators acting in L 2 (E) and Lf (E) follows from The- 
orem [3l|] and the method of continuity. For details, we refer to [2[ Theorem 4.15] 
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where it was shown that I ± EN are Fredholm operators with index zero on L 2 (E), 
(E) and 

To show that I±EN are Fredholm operators on L^(E), note that this is immediate 
from Lemma [3781 and the result for L 2 (E) when E is bounded. In case of a Lipschitz 
graph, consider the commutative diagram 



0- 



0- 



-T N 

^2 



(£)- 



I±EN 



LoCE) 



-T D 

^2 



I±EN 



-LffEl 



-0 



I±EN 



-0. 



Note that the rows are exact, i.e. the inclusion map i is injective, T is surjective 
and N(r) = R(i) = L 2 (E). It has been shown that the first two vertical maps 
are Fredholm operators. We can now apply a general technique, the five lemma, to 
deduce that / ± EN : L 2 (T?) — > L 2 (T?) is a Fredholm operator as well. For details 
concerning the five lemma for Fredholm operators we refer to Pryde [TP] , 

(2) It remains to show injectivity on L^(E). For this, it suffices to show that 
E ± L 2 H fl N ± L^ = {0} for all four intersections. If E is a Lipschitz graph, then the 
result follows from Theorem 13.31 since in this case E ± L 2 fl N ± L 2 = {0}. Assuming 
that E is a bounded Lipschitz surface, consider for example E~L 2 PI N + L 2 , and 
let F = dU = -8U, where F and U are monogenic in D~, N*(U),N*(F) G L 2 (E), 
Z7, -P — > when x —> oo, and v j F = on E. An application of Stokes' theorem 
shows that 



|F| 2 dx = / (U, v j F) da(y) - (U, 5F) dx = 0, 

and therefore F = 0. Note that the asymptotics of the Cauchy kernel E(x) shows 
that \U\ < |F| < \x\~ n and \8F\ < Ixl -1- ™ as x — > oo, so that both D~ 

integrals are convergent. 

A similar argument shows that the other three intersections also equal {0}. □ 

Corollary 3.11. If 1 < k < n — 1, then the Hardy space E ± L 2 (E; A fc ) of boundary 
traces of two-sided monogenic k-vector fields in D ± is infinite dimensional. 



Proof. Consider the Fredholm operator 



E ± : N + L 



+ T N 



±tN 



EL 



By Proposition 13.71 this maps k- vector fields to k- vector fields. Since the space of 
tangential / G N + L 2 (E; A k ) such that d%f = is infinite dimensional, the corollary 
follows. □ 



4. HlLBERT TRANSFORMS FOR fc- VECTOR FIELDS 

In this section we prove the following main result of this paper, which extends 
Theorem 12.41 to more general k- vector fields. 

Theorem 4.1. Let D C R n be a Lipschitz graph, interior or exterior domain and 
assume that 1 < k < n — 1 . 

LetU : D -> A fc R n be such that 5dU = = dSU and N t (U) , N*(dU) , N*(6U) G 
L 2 (S). If D is an exterior domain, also assume that U, dU, SU — > when x — » oo. If 
D is an interior domain, also assume that 5U = if k = 1 and dU = if k — n— 1 . 
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Then there are unique Cauchy type harmonic conjugates V\ : D — > A k 2 R n and 
V 2 : D -> A fc+2 R n suc/i f/iat 

||AW)|| 2 + ll^(^)|| 2 < ||iV*(^)ll2 + ||#. Wlh + \\N*(5U)\\ 2 . 

Before the proof of the theorem, we make some remarks. In the scalar case, we 
could apply the known results from Theorem 12.51 on L p -invertibility of the classical 
double layer potential operator to prove existence and uniqueness of Cauchy type 
conjugate functions in Theorem 12.41 For A;- vector fields, 1 < k < n — 1, Theorem 12.51 
is no longer available. In particular, we do not obtain any L p -results for p > 2, since 
the atomic estimates in the proof of Theorem 12.51 in an essential way use that the 
equation is scalar. Instead, we make use of a natural L 2 -based boundary function 
space, Lf(E), of mixed and 1 order regularity. The key observation is that both 
dll and SU are (two-sided) monogenic when U is a two-sided harmonic fc-vector 
field. Thus, in order to apply the well established L 2 -theory for BVPs, we need to 
require that N*(U), N*(dU) and N*(8U) belong to L 2 . 

Just like in the scalar case, the Cauchy type harmonic conjugate functions to a k- 
vector field U can be calculated using a generalized double layer potential operator. 
Indeed, according to (jSJ), the Cauchy integral maps 

(C±) fc - 2 Vl . 

Thus, we need to solve for h in the generalized double layer potential equation 

(tfi/Ofc = U. 

Corollary 4.2. Let D ± C R" be a Lipschitz graph, interior or exterior domain, and 
let 1 < k < n — 1. Then the range and null space of(C )k, with domain /^(E; A fc ), 
are 

/?((C±) fc ) = {U : D ± -> A fc R n ; 8dU = = d5U, N*(U), N*(5U), N*(dU) G L 2 (E)}, 
NdC^k) = {F\s ; F : D T —>■ A k R n , dF = = 8F, N*(F) G L 2 (E)}, 

wift i/ie same modifications o//?((C ± )fc) as in Theorem \4.l\ 1vhen is an exterior 
domain and when D ± is an interior domain and k = l,n — 1, and where F — > 
u>nen x — > 00 u>nen _D T is an exterior domain and F G /V((C )&). The operator 

(C*)* : (E; A fe )/ A/((C ± ) fe ) — fiftC*)*) 

zs an isomorphism. Thus, if U G RiJ^C^-)^), its Cauchy type harmonic conjugates 
are 

V 1 = (C ± ((C ± ) fe )- 1 f/) fe _ 2 ana 1 V 2 = (C ± ((C ± ) fc )- 1 f/) fc+2 . 

Proo/. If £7 G R((C ± ) fe ), i.e. if £7 = (C^/i)* for some dipole density h G (E; A fe ), 
then it follows that £7 is harmonic and N*(U), N*(6U), N*(dU) G L 2 (E). But 5dU = 
fc+ 2 = 0, so U is in fact two-sided harmonic. The converse inclusion follows 
from the existence proof of Theorem 14.11 below. 
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Clearly, if F : D T — > A fe R" satisfies dF = = 5F and iV*(F) G L 2 (S), then 
-F|e belong to the Hardy subspace for the complementary domain D T . In partic- 
ular C ± (F| 2 ) = and F| s G N((C ± ) / t). The converse inclusion follows from the 
uniqueness proof of Theorem 14.11 below. □ 

Remark 4.3. According to (JTj), the classical double layer potential is the compres- 
sion of the Cauchy integral to the subspace of scalar valued functions. In Corol- 
lary I4.2[ identifying a harmonic function with its boundary trace, we have general- 
ized this by compressing the Cauchy integral/Hardy projection E ± to the operator 
(i? acting in the subspace of k- vector fields. Other useful compressions of the 
Cauchy integral use instead the subspaces of tangential or normal multivector fields 
N ± L 2 (T,), which are relevant for BVPs. For example, consider the BVP in Exam- 
ple 13.21 consisting in finding F : D + — > AR™ satisfying the Hodge-Dirac equation, 
with a prescribed normal part g of the trace / = F\-£. Equivalently, we look for 
a Hardy function / G E + L<2 satisfying N~ f = g. Making the ansatz / = E + h, 
with h G iV~L 2 , we obtain instead a double layer type equation N~E + h = g in 
the subspace N~L 2 . As shown in ^j, the well posedness of the BVP is essentially 
equivalent to the compressed Cauchy integral N~ E + \^-l 2 being an isomorphism. 
We note that these types of compressions to N ± L 2 in general have better proper- 
ties than the compressions (-E ± )fc|L 2 (s ; A' ; ) used in Corollary 14.21 as N ± E ± \ N ± L2 are 
Fredholm operators. 

Uniqueness proof of Theorem \4 . 1\ Assume that h G L 2 (Ti; A fc ) is such that its Cauchy 
extension satisfies U = (C ± h)k = 0. We aim to prove that h G (E) and C ± h = 0, 
so that the Cauchy type harmonic conjugates V\ = (C ± ^)fc_2 and V2 = (C ± h)k+ 2 
vanish. 

Define monogenic fields 

Vi{x) := C^N-h^x) =±J E(y- x){u{y) j %)) da(y), 

V 2 {x) := C^N+h^x) =± [ E(y- x){u{y) a %)) da(y), 

for x G D ± , so that V 1 = (Vi) k -2, V 2 =jV 2 )*+2 and (Vi) k + (V 2 )k = U = by 
assumption. It follows that dV 2 = [(d + S)V 2 ) k +3 = and SV 2 = — d(V 2 )k = d(Vi) k = 
((d + 8)Vi)k+i = 0. Thus V 2 and therefore V 2 are two-sided monogenic and 

(10) T(E ± N + h) = 0. 

Similarly, it follows that T(E ± N~h) = 0. 

We first show that ( ITUl) implies that hi := N + h has regularity hi G (E). Note 
that E ± hi = ± E)hi = ~(J ± EN) hi, and consider the commutative diagram 



D 



L 2 (E)^L 2 (E), 

where the inclusion i is dense. Moreover, as explained in the proof of Theorem 13. 101 
the method of continuity shows that HEN are Fredholm operators with index zero 
on both L 2 (E) and Lf (E). Since (I ±EN)hi G (E), a general regularity theorem 
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for Fredholm operators [21 Proposition 4.16] shows that h% G Lf(E). Similarly, 
N~h G L% (E). 

We have shown that E ± (TN + h) = = E^TN-h), i.e. TN+h G E*L§ n A+Lf 
and rJV - /i G E T L^ n N~L§. Thus Theorem EHU] shows that T(N + h) = r(i\T/i) = 
0, so that ft. G (E). According to Proposition I3.7( i). we have V\ — V 2 — and 
thusC^/^0. " ' □ 

Existence proof of Theorem \4-l\ Assume that U : D ^ A fc R n is twosided harmonic 
in the sense that 8dU = = d<H7, where N*(U), N,(dU), N*(5U) G L 2 (E). If D is an 
exterior domain, also assume that U, dU, 5U — ► when x — ► 00. If £) is an interior 
domain, also assume that 5U = if k = 1 and (if/ = if /c = n — 1. 

We aim to construct V\ : D — > A fc_2 R n and V2 : D — » A fc+2 R n and a dipole density 
h G Lf (E; A k ) such that C^/i = V x + t7 + F 2 and 

||iV.(Vi)||2 + \\N*{V 2 )\\ 2 + \\h\\ D < \\N.(U)\\ 3 + \\N.(dU)\\ 2 + \\N*(8U)\\ 2 . 

(1) We first construct a tangential fc-vector field h 2 G A/" + L^(E; A fc ) such that 
dU = d(C ± h 2 ). To this end, consider the singular integral equation 

N+E^hz = N + (dU\x). 

From the assumption, dU is a monogenic field, and Proposition I3.6( i) shows that 
N + (dU\^) G L^(E). We claim that the compressed Cauchy integral 

N + E ± : N + L§(£) — > A+L^(E) 

is an isomorphism. Since it is the composition of E ± : A + L^(E) — > £ ,± L^(E) 
and A + : £ ,± L^(E) — > A + L^(E), it follows from Theorem 13.101 that it is an injec- 
tive Fredholm operator. Using the operator algebra developed for boundary value 
problems in [1], we see that 

A 2 - A(N + E + N + + N-E~N-) = (XN) 2 - (N + E) 2 

= (A + 1 + EN)N(X - 1 - EN)N, 

A 2 - 4(N + E~N + + N~E + N~) = (XN) 2 - (N - E) 2 

= (A + 1 - EN)N(X - 1 + EN)N. 

The right hand sides are seen to be Fredholm operators in L^(E) for all real A as in 
part (1) in the proof of Theorem 13. 101 Applying the left hand sides to / G A + L^(E), 
so that N~f = 0, shows that A 2 - N + E± : A+Lf (E) -> A+Lf (E) are all Fredholm 
operators. Since this operator clearly is invertible for large enough A, the method of 
continuity shows that N + E ± : A + L^(E) — ► A + L^(E) has index zero, and therefore 
is surjective, since it has been shown to be injective. 

Solving the equation, we obtain a unique h 2 G A + L^(E) such that 

N + (E±h 2 -(dU) | s ) = 0, 

where we verify that (dU)\j; G E ± L 2 V (E). It is here we need the topological assump- 
tion on E. If E is an unbounded graph, then N~L 2 fl E ± L 2 = {0}. On the other 
hand, if D + is Lipschitz diffeomorphic to the unit ball and 1 < k < n — 3, then 
(dU)\^ G E ± L^(H) according to Lemma ESI lik — n — 1 then dU = G E ± L 2 H (Y) 
by assumption. Finally, if k = n — 2 then Stokes' theorem shows that 

/ v a (d[/)dcr = ± / d(dU)dx = 
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since d 2 = 0. (Note that Stokes' theorem is applicable in D since \dU\ < \x\ n as 
x -> oo.) Thus E ± h 2 -{dU)\ 11 G E ± L^nN~L^, which equals {0} by Theorem EH 
In either case, we conclude that 

dU = d(C ± /i 2 ), for some h 2 G iV + Lf(E), 

such that h 2 = r/i 2 and \\N,{U)\\ 2 + \\N*(dU)\\ 2 > \\h 2 \\ R « ||/i 2 || 2 + ||/i 2 || 2 . 
A similar argument proves that 

5U = 5(C ± /n), for some h x G N~L^{T), 

such that 1 1 /ix 1 1 2 + \\Th t \\ 2 < \\N*{U)\\ 2 + \\N*(8U)\\ 2 . 

(2) To construct Cauchy type harmonic conjugates V\ and V 2 , write V, := C^hi, 
i = 1,2. Then 

c?(C/ - (Vi)* - (V 2 )fc) = <«/ - - dtf = 0, 
S(U - (yi) k - (V 2 ) k ) = 8U-6U-0 = 0, 

so that (U — (Vi)k — (t^)fe)|s G E ± L 2 . Thus, defining a dipole density 
h:=(U- (t4) fc - (y 2 ) fc )| E + /ix + /i 2 G A fe ), 

Vi := (Vi)fe_2 and V 2 := (V 2 )fc+2 gives 

C ±h = U- (Vi)* - (V 2 ) fc + Vi + V 2 = 14 + [/ + V 2 . 

Since ||A^i)ll2 + ||AUy 2 )|| 2 < ||JV,(E/)|| 2 + |Ml2+IIM2 < ||iV,(Z7)|| 2 +||JV,(d^)|| 2 + 
||JV.0H7)|| a and ||/i|| D < || (C/ - {V x ) k - (V 2 ) fc )| S || 2 + ||/n|| D + ||/i 2 || c < ||iV*([/)|| 2 + 
||iV*(e#7)||2 + \\N*(5U)\\ 2 , the proof of Theorem O is complete. □ 

5. Other types of harmonic conjugates 

Recall from the discussion in the introduction that if U : D — > A fc R n is a two-sided 
harmonic fc-vector field, two fields V\ : D — »■ A fc ~ 2 R™ and V 2 : Z) — > A fc+2 R n are said 
to be conjugate to U if (d + S)(Vi + U + V 2 ) = 0. As noted, further conditions 
need to be imposed on Vi, V 2 for this problem to be well-posed, i.e. for Vi to be 
(essentially) unique. Theorems 12.41 and 14.11 show that the problem becomes well 
posed if Vi are required to be Cauchy type harmonic conjugates. The following 
proposition expresses this condition as a boundary condition. 

Proposition 5.1. Assume that V\ and V 2 are harmonic conjugates to U in D = D + , 
with N*(U), N*(Vi), N^(V 2 ) G L 2 (E). Then they are of Cauchy type if and only if 
there exists U~ : D~ — > A fc R™ with harmonic conjugates V{~ : D~ — > A fc ~ 2 R™ and 
V 2 : D~ -> A fc+2 R" in D-, with N,(U~), iV*(Vf ), N*(V 2 ) G L 2 (E) (and decay at 
infinity when D~ is an exterior domain), such that 

Vi| s + Vf| s = and V 2 \ s + V 2 ~| s = 0. 

Proof. Recall that by definition, Vi are of Cauchy type if there exists h : S — ► A fc R™ 
such that C + /i = V± + U + V 2 . In this case, Vf~ + C/~ + V 2 ~ := C~/i has the required 
properties since E + h + = h. Conversely, if V{~ + U~ + V 2 has the required 
properties, let h := C/| s + C/-| S = (Vi + U + V 2 )| E + (Vf + U~ + V 2 -)| E : £ -> A fc R n . 
Then C+/i = V x + £/ + V 2 . □ 
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Note that since Cauchy type conjugates are defined in terms of the Cauchy in- 
tegral, which concerns the interplay between monogenic fields in D + and D~~, the 
boundary condition above is a transmission problem, i.e. a jump relation between 
pairs of monogenic fields in D^. 

We end this paper with a construction of harmonic conjugates which differ from 
the Cauchy type ones in general. To avoid technicalities, we shall only consider 
interior domains, i.e. D is assumed to be Lipschitz diffeomorphic to the unit ball. 
(Unlike the situation for the Cauchy type conjugates, this second construction does 
not involve the complementary domain D~ .) 

Consider the exterior and interior derivative operators d and 5 in ([2]). These are 
formally (anti-) adjoint in the sense that j Rn (8F(x), G(x))dx = — J Rn (F(x), dG(x))dx 
for all F,G G C^°(R n ; A), which is a consequence of the definition (jlj) of the interior 
product. On the domain D, on the other hand, an application of Stokes' theorem 
shows that 

/ (5F,G)dx+ [ (F,dG)dx= [ (vjF,G)da= [ (F, v a G)do~. 
Jd Jd Jt, Jt, 

Thus, in order to make d and —5 to be adjoint operators in L 2 (D; A), one needs to 
impose either tangential boundary conditions on 5 (so that v j F\j^ = 0) or normal 
boundary conditions on d (so that u a G\s = 0) . If these boundary conditions are 
imposed on the operators in a suitable weak sense, then one obtains two pairs of 
densely defined, closed and adjoint operators (—8, d) and (— 8, d) in L 2 (D; A), where 
the domains of the operators are 

D(8) = {F e L 2 ; 8F G L 2 }, D{8) = {F G L 2 ; 8F G L 2 ,u j F| s = 0}, 
D(d) = {GeL 2 ; dGe L 2 }, D(d) = {GeL 2 ■ dGeL 2 ,uA G| s = 0}, 

i.e. 8 is 5 with tangential boundary conditions and d is d with normal boundary 
conditions. For further details, we refer to Section 4]. The following proposition 
summarizes relevant facts from the theory of Hodge decompositions that we need 
below to construct harmonic conjugates. 

Proposition 5.2. Assume that D C R™ is Lipschitz diffeomorphic to the unit ball. 

(i) If F G L 2 (D; A k ) satisfies dF = 0, and 2 < k < n, then there exists a unique 
V G L 2 (D] A*" 1 ) such that 8V = 0, v j V| E = and 

F = dV. 

If k = 1, then such V exists, but is unique only up to constant scalar func- 
tions. 

(ii) If F G L 2 (D; A fc ) satisfies 8F = 0, and < k < n — 2, then there exists a 
unique V G L 2 (D; A k+1 ) such that dV = 0, v a V\y, = and 

F = 8V. 

If k = n — 1, then such V exists, but is unique only up to constant n-vector 
fields. 

In terms of operators, this result means that 

d-.N&A*- 1 ) — >N{d;A k ), and 5 : N(d; A k+1 ) — >N(5;A fc ), 
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are injective (except when k = l,n — 1 respectively) and surjective unbounded 
operators. This follows from [51 Theorem 1.3] together with Poincare's lemma. Here 
N(T; A fc ) := N(T) fl L 2 (D; A k ). 

Definition 5.3. Given a two-sided harmonic field U : D — > A fc R n , i.e. <Jd?7 = = 
c/<5[7, we say that V\ : D — > A fc ~ 2 R n and V2 : -D — > A fc+2 R n are Hodge type harmonic 
conjugates to U if (cf + 5) (Vi + £/ + V2) = and if f j Vi\z — and v a V^Ie = 0. 

Our main result in this section is the following theorem, which shows that the 
boundary condition imposed on Hodge type harmonic conjugates yields a well posed 
problem. 

Theorem 5.4. Assume that D C R n is Lipschitz diffeomorphic to the unit ball, let 
< k < n and let U e L 2 (D; A k ) be such that dU, 5U e L 2 (D; A) and 5dU = = 
dSU. If k = n — 1, assume that dU = 0, and if k = 1 assume that 5U = 0. Then 
there exists Hodge type harmonic conjugates Vi,V 2 to U in D such that ||Vi||z, a (.D) ^5 
||<5£^||l 2 (d) and \\V 2 \\l 2 (d) % ||^||l 2 (-D)- The conjugates are unique, except if k = 2, 
when Vi is unique modulo constants, and if k = n — 2, when V 2 is unique modulo 
constants. 

Proof. Apply Proposition I5.2( i) to F — SU, and Proposition I5.2( ii) to F — dU. □ 

Remark 5.5. (1) In the complex plane, when n = 2 and k = 0, harmonic con- 
jugates in general are unique modulo constants. In particular Hodge type con- 
jugates coincide with Cauchy type conjugates, modulo constants. Note that all 
V — V 2 : D — > A 2 R 2 are normal on the boundary, since v a V|s = 0. 

(2) For general domains D, Cauchy type harmonic conjugates and Hodge type 
conjugates will not coincide in general when n > 3, not even when k = 0. To see 
this, note that there is no reason for the Cauchy type conjugate 

V{x) = [ (E(y - x) a v(y))h{y)do-(y), h : E — > A°R" 
Jt. 

to satisfy v a V\s = 0. However, they do coincide, for all < k < n, when D is a 
sphere. To see this for the unit sphere, note that the normal vector u(y) in this case 
is y, so 

y a h(y) da(y) 
o- n _i|y - x\ n 

is normal on E. Similarly, the Cauchy type conjugate V\ is seen to be tangential 
on E in this case. Hence they coincide with the Hodge conjugates by uniqueness in 
Theorem 15.41 (modulo constants when k = 2, n — 2). 
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